First Study of the Negative Binomial Distribution Applied to Higher Moments of 
Net-charge and Net-proton Multiplicity Distributions 
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A study of the first four moments (mean, variance, skewness, and kurtosis) and their products (kct 2 
and Sa) of the net-charge and net-proton distributions in Au+Au collisions at ^/snn = 7.7-200 GeV 
from HIJING simulations has been carried out. It is seen that a Poisson does not effectively describe 
the actual distributions of positive and negative particles (or protons and anti- protons). A discrete 
probability distribution that effectively describes the raw distributions is the negative binomial (or 
binomial) distribution (NBD/BD). The NBD/BD have been used to characterize particle production 
in high-energy particle and nuclear physics. Differences between ko 2 and the Poisson assumption of 
a factor of four (for net-charge) and 40% (for net-protons) can be accounted for by the NBD/BD. 
This is the first application of the NBD/BD to describe the behavior of the higher moments of 
net-charge and net-proton distributions in nucleus-nucleus collisions. 

PACS numbers: 25.75.-q, 25.75.Gz 



INTRODUCTION 

Colliding heavy ions is an essential experimental tool 
for exploring the nature of the deconfinement phase tran- 
sition. The phase transition line between hadronic and 
partonic matter in the phase diagram of quantum chro- 
modynamics (QCD) is frequently drawn as a function of 
baryon-chemical potential (hb) and temperature. Heavy 
ion collisions at a particular center-of-mass (CMS) en- 
ergy essentially probe a singular value (small range) of 
fiB- Baryon-chemical potential can be changed by ad- 
justing the CMS collision energy. The Relativistic Heavy 
Ion Collider (RHIC) recently embarked on a program to 
study the QCD phase diagram by colliding Au-nuclei at 
CMS collision energies (^/snn) of ^/snn = 7.7, 11.5, 19.6, 
27, 39, 62.4, and 200 GeV Q. This 'energy scan' provides 
information at different values of {1b from ss 20-400 MeV 
Lattice QCD (LQCD) and experimental data sug- 
gests that if there is a phase transition at high tempera- 
tures close to ^ B = 0, it appears to be a smooth crossover 
Q. Additional predictions from LQCD (dependent on 
the number of quark flavors and their masses) indicate 
the possible existence of a first order phase transition at 
larger values of /ib- If this is the case, where the first or- 
der phase transition line and crossover phase transition 
region meet could be a critical point [1] . 

From previous studies of the thermodynamics of phase 
transitions, it has been seen that a system undergoing 
a phase transition that passes through a critical point 
can exhibit enhanced fluctuations (e.g. critical opales- 
cence) Lattice QCD studies have shown that 
near the phase transition the fourth-order quark num- 
ber susceptibilities of conserved quantities (baryon quan- 
tum number, strangeness quantum number, and charge) 
change rapidly [7] . Because experiments cannot measure 
all produced baryons, strangeness, or charge, the higher- 



order moments of the following proxy distributions are 
measured for each conserved quantity (baryon quantum 
number, strangeness quantum number, and charge): net- 
proton, net-kaon, and net-charge, respectively. In this 
case, the first through fourth-order moments correspond 
to the mean, variance, skewness, and kurtosis, respec- 
tively. There have been estimates that the skewness (S) 
and kurtosis (k) should be sensitive to critical fluctu- 
ations generated by a phase transition trajectory that 
passes through (or close to) a critical point. The higher 
moments are proportional to powers of the correlation 
length (£) with the skewness (third moment) propor- 
tional to £ 4,5 and the kurtosis (fourth moment) propor- 
tional to £ 7 0. 

The moments of the (e.g. net-charge) distribution con- 
tain a trivial system volume dependence. Constructing 
products and ratios of these moments can cancel the 
intrinsic volume dependence. Two examples of these 
volume-independent products are Sa and na 2 . These 
ratios of moments can also be expressed in terms of cu- 
mulants: Sa = ^ and na 2 = Quantitatively, their 
sensitivities to the correlation length (£) are Sa oc £ 25 
and no- 2 cx £ 5 ||. 

Several measurements of the higher-moments of net- 
proton, net-kaon, and net-charge distributions have been 
carried out Tol- 12\. These measurements are frequently 
compared to a baseline that assumes Poisson statistics. 
Deviations from this Poisson baseline have been ob- 
served. However, the deviations from the Poisson base- 
line are not as large as those expected from strong en- 
hancement due to the presence of a QCD critical point. 
One estimate of na 2 for the net-proton distribution from 
a QCD-based model that includes a critical point is 
kit 2 « 2.5, 35, or 3700 for a critical point at ^/snn = 200, 
62.4, and 19.6 GeV, respectively [10| . Other predictions 
concluded that the sign of the skewness and/or kurtosis 
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will change at a critical point 13|, ll4[ . This behavior has 
been observed from lattice QCD studies that have a crit- 
ical point, whose location in energy and /iB is dependent 
on parameters such as the critical temperature (T c ) and 
the lattice spacing used [l5j . 

Before interpreting the results from the higher- 
moments studies in terms of critical point phenomena, 
the contribution from known physics must be quanti- 
fied. It is well known that experimental multiplicity 
distributions are better described by the (negative) bi- 
nomial distribution (NBD/BD) [T^ED]. The NBD/BD 
and Poisson distributions are all discrete probability dis- 
tributions. The Poisson distribution is a limiting case 
of the NBD/BD where the mean and the variance are 
equal. For a NBD/BD distribution, the variance is 
larger/smaller than the mean. There are several sources 
of correlation in particle production that can cause devi- 
ations from pure Poisson statistics: e.g. charge conserva- 
tion, correlated production of particles from one source 
(e.g. strings or clusters), etc. Therefore, it is impor- 
tant to verify the experimental results compared to not 
only a Poisson distribution, but to other valid baselines 
that better describe the multiplicity distributions. The 
NBD/BD baseline is one of these cases. 



and weighted by the number of events in each bin to 
remove any centrality bin-width effects 24 , 2{| . 



RESULTS AND DISCUSSION 

Figure [T] shows the raw histogram for positively 
charged particles from ^/snn = 200 GeV Au+Au col- 
lisions (left) and protons from ^/snn = 7.7 Gev Au+Au 
collisions (right) in a representative centrality bin (10- 
12.5%) from HIJING simulations. Two baseline distri- 
butions are calculated from the information in the his- 
tograms: the solid lines are a negative binomial distribu- 
tion (NBD) and the dashed lines are a Poisson distribu- 
tion. The only input to the Poisson curve is the mean of 
the number distribution, while for the NBD the mean and 
variance of the number distribution are required. For the 
positively charged particles, the x 2 /NDF for the Poisson 
is 91900/162 («567) and for the NBD is 195/162 (1.2). 
For the protons, the \ 2 /NDF for the Poisson distribu- 
tion is 2577/38 (68), while for the NBD the \ 2 /NDF is 
34.6/38 (0.91). In both cases, the NBD is a better de- 
scription of the data than the Poisson distribution. The 
same is true for negatively charged particles and anti- 
protons. 



ANALYSIS 

For this study, minimum bias Au+Au events at ^/snn 
= 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV were created 
with the default HIJING (vl.383) event generator [22 1. 
Distributions of positive and negative particles and pro- 
tons and anti-protons were produced. From each distri- 
bution, the mean value and the variance were extracted. 
From this information, the higher order cumulants for 
each individual distribution could be calculated. From 
these, the combined cumulants of the net-charge and net- 
proton distributions were calculated. The additive prop- 
erties of cumulants means that even order cumulants are 
added and odd order cumulants are subtracted (also for 
the NBD [23[) to calculate the required cumulant for the 
net-charge/proton distribution. 

The centrality is defined using charged particle mul- 
tiplicity in the region 0.5 < \r}\ < 1.0. The measured 
charged particles and identified protons/anti-protons are 
restricted to the pseudorapidity range |?7| < 0.5. For 
this measurement, this is required to prevent auto- 
correlations between the particles used for the measure- 
ment of the moments of the multiplicity distributions and 
those used for centrality determination. This method 
has been used in other analyses of particle correlations 
in pseudorapidity [24|]. For net-charge studies, the trans- 
verse momentum range for particles in this analysis was 
0.2 < pt < 2.0 GeV/c and for net-proton studies the 
range px > 0.4 and p < 3.0 GeV/c was used. 

The moments are calculated in each multiplicity bin 
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FIG. 1. Raw distributions (filled circles) of positively charged 
particles from ^/snn = 200 GeV Au+Au collisions (left) and 
protons from v / snn = 7.7 Gev Au+Au collisions (right) in the 
10-12.5% most central collisions from HIJING. Two baseline 
distributions are calculated using the mean value (for Poisson) 
(solid line) or the mean and variance (for negative binomial 
distribution (NBD)) (dashed line) of the raw distributions. 
The x 2 /NDF of the Poisson and NBD compared to the raw 
distributions are also included. 

The centrality dependence of kct 2 (C4/C2) for the net- 
charge and net-proton distributions calculated from sim- 
ulated Au+Au collisions at ^/snn = 7.7, 19.6, 39, and 200 
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GeV from HIJING are shown in Figures 2(a) 



The two bins with largest N pavt correspond to 0-5 and 
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5-10% most central collisions, while the other seven bins 
are 10% wide bins incorporating the remaining 10-80% ol 
the hadronic cross-section. The results from HIJING for 



net-charge (Figure 2(a) I demonstrate a centrality inde- 
pendent trend for na 2 , with similar average magnitudes 
at all energies from ^/snn = 7.7-200 GeV of between w2- 
2.5. The Poisson values for na 2 are always exactly equal 
to one. For the NBD baseline the values are also inde- 
pendent of centrality, except for the most central bin, 
which shows an anomalous decrease. This decrease in 
the most central bin is also observed in the net-proton 



results from HIJING (Figure 2(b)) and in the data anal- 
ysis (loj . This downward fluctuation in the most central 
bin seems to persist in both hadronic models and data, 
which indicates that it is almost certainly a consequence 
of the analysis method. This effect is not seen in the 5- 
10% bin, which generally follows the overall trend of the 
other centralities. 

The NBD baseline is also larger than the measured 
na 2 from HIJING. One possible explanation for this dif- 
ference is that particle production in HIJING does not 
always follow a NBD, but contains additional contribu- 
tions. Though the NBD is larger in magnitude, the qual- 
itative trend of na 2 for net-charge is reproduced 

Figure 2(b) shows that no 2 (C4/C2) for the net-proton 



distribution from HIJING is also approximately central- 
ity independent and have values between sal. 0-1. 5. As 
before, the Poisson baseline for no 2 is always exactly 
equal to one. The NBD baseline provides an excellent 
description of the centrality dependence of na 2 for net- 
protons. 

In all cases, the statistical error bars are likely to under 
estimate the true magnitude of the error. The statistical 
errors on the variables Sa and na 2 have a non-trivial 
dependence on the shape of the underlying distribution. 

The statistical error is dependent on the variance of 
these underlying distributions (26|. Therefore, for the 
same number of events when the number of particles in 
the analysis is decreased (e.g. net-protons instead of net- 
charge) the size of the statistical error also decreases be- 
cause the variance of the net-proton distribution is nar- 
rower than that of the net-charge distribution. As such, 
this method cannot be a wholly unbiased estimation of 
the statistical error, particularly for cases where the vari- 
ance of the underlying distribution is small solely due to 
limited numbers of accepted particles. 

Figures [3] and H] shows the incident energy dependence 
of Sa (C3/C2) and kg 2 {C 4 /C 2 ) for HIJING simulated 
central 0-5% Au+Au collisions at ^/snn = 7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV for net-charge and net- 
proton distributions, respectively, na 2 for both the net- 
charge and net-proton distributions exhibit no overall en- 
ergy dependence. Sa demonstrates a smooth decrease 
from y/sNN = 7.7-200 GeV. For net-charge (net-proton) 
measurements, Sa is intrinsically related to the nega- 
tive/positive (anti-proton/proton) ratio. Also plotted in 
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FIG. 2. Centrality depen dence of kct 2 (C4 /C2) ( open circles) 
for the net-charge (2(a) I and net-proton ( 2(b) l distributions 
from HIJING simulated Au+Au collisions at ^/snn = 7.7, 
19.6, 39, and 200 GeV. Also shown are predictions assum- 
ing the underlying distribution is Poisson (dashed line) or a 
negative binomial distribution (NBD) (solid line). 



Figures [3] and |4] are the Poisson (dashed line) and NBD 
(solid like) baselines for both Sa and na 2 . It is clear 
from the comparison that the NBD provides a better de- 
scription of the energy dependence of Sa and na 2 than 
the Poisson baseline for both net-charge and net-proton 
distributions. Some of the values of na 2 for 0-5% most 
central HIJING events do fall below the NBD baseline. 
However, as discussed above this is a systematic artifact 
of the analysis method that occurs in the 0-5% bin. Com- 
paring the 0-5% and 5-10% bin (which should nominally 
contain the same physics) in Figures 2(a) and 2(b) it is 



shown that the 5-10% bin does not exhibit these strong 
excursions from the overall centrality trend. 



4 




FIG. 3. Energy dependence of kct 2 (C4/C2) and Sa (C3/C2) 
for net-charge distributions from HIJING simulated central 
Au+Au collisions (0-5%). Also shown are predictions assum- 
ing the underlying distribution is Poisson (dashed line) or a 
negative binomial distribution (NBD) (solid line). 
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FIG. 4. Energy dependence of na 2 (C4/C2) and Sa (C3/C2) 
for net-proton distributions from HIJING simulated central 
Au+Au collisions (0-5%). Also shown are predictions assum- 
ing the underlying distribution is Poisson (dashed line) or a 
negative binomial distribution (NBD) (solid line). 

SUMMARY AND CONCLUSION 

The centrality and energy dependence of the moments 
(mean, variance, skewness, and kurtosis) and their prod- 
ucts (k<t 2 and Sa) of net-charge and net-proton mul- 
tiplicity distributions in Au+Au collisions simulated by 
HIJING at ^/snn = 7.7, 11.5, 19.6, 27, 39, 62.4, and 
200 GeV have been studied. These results have been 
compared to two cases where the underlying distribution 
is assumed to be either a Poisson (uncorrelated parti- 
cle production) or a negative binomial/binomial distri- 
bution (NBD/BD). The negative binomial distribution 
has been shown to accurately reproduce features of mul- 



tiplicity distributions in high-energy elementary (pp) and 
nucleus-nucleus collisions. It has been demonstrated that 
the Poisson is a poor approximation of the raw distribu- 
tion for both positive (and negative) particles or protons 
(and anti-protons). For the first time it is shown that 
the NBD/BD is a better approximation of the raw dis- 
tributions and provides a comprehensive description for 
the main centrality and energy dependence trends of the 
moments products na 2 and Sa. This unique insight in- 
dicates that relatively small deviations from the Poisson 
assumption, which are well described by the NBD/BD 
assumption, cannot provide an unambiguous indication 
that Ka 2 and Sa are sensitive to the quark/gluon-hadron 
phase transition or any potential QCD critical point. 
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